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Abstract 

Proximity effect systems in superconducting films can be modeled by a one-dimensional 
Schrodinger equation. Several systems are studied using Dirichlet and Neumann boundary con- 
ditions. It is observed that the two boundary conditions have a dramatic effect on the lowest 
eigenstate allowed in these systems, and points to unusual behavior for solutions of Schrodinger's 
equation in certain potential wells and proximity effect systems. 



I. INTRODUCTION 

The superconducting proximity effect is the change in the superconducting transition 
temperature of a superconducting film when another metallic film with either a lower or 
zero superconducting transition temperature is deposited onto it, resulting in a proximity 
effect bilayer. There is much interest in these systems as transition edge calorimeters, such 
as the Mo/Cu bilayers being developed for detecting radiation.^ Predicting the transition 
temperature of such systems is usually done with the de Gennes-Werthamer theory^ or 
the theory of UsadelP These theories are difficult to introduce at the level of introductory 
quantum mechanics or solid state physics. However, it is possible to make a one-to-one 
correspondence between the problem of proximity effect in layered films and the quantum 
mechanics textbook example of the energy levels and wave functions of a particle in a one- 
dimensional potential wellP 

Reference |3| derives the standard equation of proximity effect systems in the dirty limit 
(electron mean free path much less than the superconducting coherence length) given by (in 
one dimension) 



where x( z ) = VKV^ + z /2) ~~ "0(1/2), with ip(x) here being the diGamma function, 5 A(x) 



superconductor or normal metal at point x, 9d and T c (x) are the Debye temperature and 
the superconducting transition temperature for the material of interest at x in isolation, 
and T c is the superconducting transition temperature of the composite structure. Typically 
T c (x) is assumed to be uniform within a particular layer. This equation is the differential 
equivalent of an integral equation for the proximity effect, and is usually much easier to 
solve. Typically it is solved by expressing the self energy as a Fourier series and looking 
for the allowed wavevectors. As shown in Reference G2, the above equation will predict the 
superconducting transition temperature of a proximity effect bilayer, however the solutions 
of this equation are difficult to present in undergraduate courses. 

Reference El discussed how equation Q looks like the one- dimensional time independent 
Schrodinger equation, 
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where ip(x) in this equation (and in the rest of the paper) is the wavefunction which is 
related to A(x), h 2 /2m is related to £ 2 , V(x) is related to hi{9 D /T c {x)), and E is related to 
ln(#£i/T c ). (The simplifying assumption is made that the quantity £ is the same for both 
layers, which in the Schrodinger equation is equivalent to assuming that the mass of the 
particle is the same in both layers.) Because solutions for constant potentials are easy to 
obtain with Eq. Q, it can be used to help students see how the transition temperature 
of proximity effect systems changes as layer thicknesses and potentials are modified. The 
purpose of this paper is to give some example solutions for various geometries and different 
boundary conditions and discuss how the solutions differ from what is typically expected. 

II. DESCRIPTION OF THE SPECIFIC PROBLEM 

Although solutions to Eq. ^ will not give precise answers to the values of T c in a real 
system of metal layers, they will reproduce the overall behavior, provided we recognize 
several limitations. First, the energy eigenvalues of Eq. ^ are inversely related to the 
transition temperatures. The lowest energy eigenstate corresponds to the highest transition 
temperature. Because ip is related to the function A in Eq. Q, we also need to consider the 
boundary conditions. Although A and its derivative are not continuous,^ we will use the 
continuity of ip and dip/dx to solve a simpler problem, and see that it still reproduces almost 
the same equation used in Ref. SI (The actual boundary conditions used in solutions of the 
superconducting proximity effect are given in Ref. El) Another issue is what to do with the 
end surfaces, which in a proximity system are the metal/vacuum interfaces. In the usual 
particle in the box problem we exclude the possibility of the particle being outside the box, 
and we force ip to go to zero at a boundary. However, for A in layered superconductors, 
the correct boundary condition at a metal/vacuum boundaries is dA/dx = 0, which implies 
that the Cooper pairs are reflected off these boundaries back into the metal layer™ We will 
compare similar geometries with two boundary conditions at the extreme boundaries, a 
Dirichlet type (if) = 0) and a Neumann type (dip/dx = 0). 

To mimic a bilayer system the potential is as shown in Fig. fit a), where the zero potential 
for < x < d represents a material with a high transition temperature, and the finite value 
of V for d < x < 2d represents a material with a lower transition temperature. (We can set 
the lower potential to any value we choose to adjust the zero of the energy.) Because the 
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probability density is assumed to vanish outside these two layers, there are infinite barriers 
at the edges. Any sequence of high and low transition temperature layers can be mimicked, 
as seen in Fig. Rib) for a 2 period proximity system, and it would be easy to vary the ratio of 
the thicknesses of the high and low potential. For the purpose of this paper equal thickness 
layers will be used. 

Rather than changing the thickness of the layers, the potential V will be varied and 
eigenvalues and eigenstates found. There is a one-to-one correspondence between the two, 
for in our solutions all energies and potentials will be given in units of h 2 /(2md 2 ). For 
example for a bilayer with width 2d, the lowest energy eigenstate would be 7r 2 /4 with V = 0, 
the lowest one for an infinite square well of width 2d (with a Dirichlet boundary condition). 
So changing V is equivalent to changing d and vice versa. 

We have to interpret the meaning of ip for the case of the proximity effect. As stated in 
Ref. [2]when A(r) depends on only one space coordinate, x, the energy gap E is equal to the 
minimum value of |A(x)| in the sample. This role of A means that it can not be allowed to 
go to zero or be negative,^ and hence if) must be positive definite. In addition, the transition 
temperature of a bilayer system will always be between the transition temperatures of the 
two materials so that only solutions with energy eigenvalues satisfying < E < V will be 
studied. There are valid solutions with E > V, but they will not shed light on the proximity 
effect. 



III. ONE PERIOD 

We first look at a simple bilayer, as seen in Fig. [l^a), which in the absence of the potential 
V would have the lowest eigenstate E = tt 2 /A. We seek a solution that satisfies < E < V. 
Because the potentials are constant and < E < V, we know if) can be written in the 
low potential layer (layer 1 with the higher transition temperature) as if>i(x) = aci cos(kx) + 
Pi sm(kx) where k = y 2mE/h 2 ; in the high potential layer (layer 2 with the lower transition 

temperature) ip2(x) = «2 cosh(ga;) +P2 sinh(gx) where q = J 2m(V — E)/Ti 2 . The continuity 
of if) and dijj/dx will be enforced at x = d, but we first have to consider the boundaries at 
the extreme edges, which represent the metal/vacuum interface. 

We first look at the Dirichlet boundary condition, which forces if) to be zero at x = and 
x = 2d. For this case the first layer will have ipi(x) = Pisin(kx) and for the second layer 
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^2 (#) = 02 sinh[g(2<i — x)]. Continuity of ip at x = d gives /3isin(A;d) = sinh(gd), and 
continuity of the derivative of ■?/> gives kf3\ cos(kd) = —qfo cosh(gd). The ratio of these two 
conditions gives, 

qtan(kd) + ktanh.(qd) = 0. (3) 

The roots of Eq. (|3| and all others were found graphically^ recall that the values of E and 

V are given in units of h 2 /2md 2 . For this case there are no eigenstates with energies less 
than V until V ^ 4.12; the energy eigenvalues are shown in Fig. [2j There is a lower energy 
eigenstate for V < 4.12, but it is greater than V. (For V < 4.12 the lowest solutions are 
given by solving qtan.{kd) + ktan(qd) = 0, whose solutions are plotted in Fig. |2j) The limit 

V — > oo is equivalent to an infinite well with width 1/2, or E = tt 2 « 9.870. 

For the Neumann condition the wavefunction must have a zero slope at x = and 
2d. The form that satisfies these conditions is ipi(x) = ot\ cos(kx) in the first layer and 
^(x) = 0.2 cosh(g(2<i — x)) in the second layer. The continuity of ip and its derivative gives 
«i cos{kd) = ot 2 cosh(grf) and — ka\ sm{kd) = — qa 2 sinh(g<i). The ratio gives the energy 
eigenvalue equation, 

ktem(kd) = gtanh(grf), (4) 

which is the same as Eq. (14) in Ref. [3] and similar to Eq. (4.16) in Ref. El In contrast to 
the eigenvalues for the Dirichlet boundary condition, the energy eigenvalues satisfy E < V 
for all values of V. If V — 0, which would represent a layer of material with the higher 
transition temperature with width 2d, T c of the bilayer would be the same as that of the 
higher transition temperature layer, which we represent by E = 0. The energy eigenvalues 
are shown in Fig.[2]in comparison to those for the Dirichlet condition. We see, consistent with 
the exact proximity effect equations, that as d increases (or V decreases), E approaches zero, 
and the predicted transition temperature approach the higher T c material; as d decreases, 
the predicted transition temperature decreases. 

The corresponding wavef unctions are compared in Fig. [3] for V — 5. As expected the 
Dirichlet condition puts the bulk of the probability in the lower well. In the Neumann 
case the shape is what is expected for a standard bilayer proximity effect system, with the 
superconducting order parameter decaying into the non-superconducting system, as shown 
in Ref. [2j The solutions in Fig. [3] have been properly normalized, but for the proximity effect 
(Neumann) case there is no clear normalization condition, because the minimum of the gap 
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function would be the measured energy gap in the systemP 



IV. TWO PERIODS 

Next we turn to the effect multiple layers will have on the lowest eigenvalue of our 
quantum potential, or the transition temperature of a multiple layered structure. Consider 
the potential is shown in Fig. [iXb), which would represent a two period bilayer. For Dirichlet 
boundary conditions if) must be chosen in the first layer as ipi(x) = (3\ sin(fcx) and in the last 
layer as ip^(x) = 04 sinh[g(4x — d)]. For the interior layers (2 and 3) it is convenient to chose 
^2 (#) — a 2 cosh[g(2<i — x)] + /3 2 sinh[g(2<i — x)] and ip 3 {x) = a 3 cos[A;(x — 2d)] + /3 3 sin[k(x — 
2d)]. Continuity of ip and its derivative at the three interior boundaries gives six relations, 
which can be reduced to eliminate the coefficients in tp giving the equation for the energy 
eigenvalues, 

[2kq + (q 2 - k 2 ) tan(fcd) tnnh(qd)]u D = 0, (5) 

where ud = [qta,n(kd) + fetanh(gd)] is the solution of the one period Dirichlet condition 
problem in Eq. Thus, solutions to the one period case are also solutions to the two 
period case. 

For the moment we will assume that ub 7^ 0, so that the two period eigenvalue equation 
becomes 

2kq + (q 2 - k 2 ) tan(Jfed) tanh(gc/) = 0. (6) 

The two period energies given by removing the one period solutions are shown in Fig. [2] and 
are lower than those of the one period Dirichlet solution. In Fig. [4] the two wavefunctions 
for the two period case are shown, one using the energy eigenvalue of the one period case, 
Eq. ([3]), and the other the energy eigenvalue for the two period case with the one period 
solution removed, Eq. ([6]). In both cases V — 5. The bulk of the probability for the lower 
energy state is in the third layer. For the higher energy state the bulk is in the first layer, 
where the Dirichlet condition in combination with the one period solution (which creates a 
node at x/d = 2) forces a smaller wavelength, resulting in a larger energy. The result is what 
we would expect, namely, if the system size increases, the energy eigenvalues of the system 
decrease. The interesting feature is that there is a lower energy eigenvalue that appears 
when the system goes to two periods, but the original energy eigenvalue does not change as 
the system size increases. 
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For the Neumann condition the only modification for two periods is that the values of 
ipi(x) and ipi(x) are changed to a\ cos(fcr) and a 4 cosh[g(4x — d)], respectively, to enforce 
zero slope for ip at the edges. This change and continuity of if} and its derivative in the 
interior three boundaries give the energy eigenvalue equation, 

[2kq + (q 2 - k 2 ) ttm(kd) tan\i(qd)]u N = 0, (7) 

where un = [gtanh(grf) — ktan(kd)], which is the same as Eq. Q. Equation is almost 
the same as derived in Ref. [7] for the full DeGennes-Werthamer solution for a two period 
bilayer; the interest here is how the eigenvalues and eigenfunctions compare. As was seen for 
Dirichlet boundary conditions, the one period solution is also a solution to the two period 
case, so the repetition of the single period solution is not due to the boundary condition. As 
before, can be factored out to give the reduced eigenvalue equation which is the same as 
that for the Dirichlet boundary condition, Eq. ([6]). Comparison of these energy values to the 
one period Neumann condition is also shown in Fig. [2j but here the energies for Eq. Q are 
higher than those for the one period Neumann condition, Eq. Q. Thus, in the Neumann 
case, increasing the size of the system does not lower the lowest energy eigenvalue. As can 
be seen m Fig. § the one period solution enforces a zero slope at the x = 2d boundary, and 
creates a reduced version of the first half for the second half. The other solution has a node 
in the wavefunction which is not allowed for A(x) for the proximity effect. 

For the two period case there is an allowed solution which repeats the one period solution, 
but whose energy compares with the lowest eigenvalue very differently. For the Dirichlet 
condition it is higher, and for the Neumann condition it is lower than the common energy 
eigenstate, whose energy does not depend on the boundary conditions. Thus, although 
there are energy eigenstates that are similar for the two boundary conditions, there are also 
differences that can be used to show the importance of boundary effects. 

V. THREE PERIODS 

In this case the potential has three repeating sections going from x/d = to 6, and 
the complexity of the calculations increases substantially. We choose Dirichlet conditions 
on the ends and label the layers as 1,. . . , 6. We must choose ipi{x) as before, but choose 
iPq(x) = p 6 smh[q(6d — x)]. We use the same values of fafo) an d ip3{x), with ipi(x) = 
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«4 cosh[g(4c/ — x)] + fl± sinh[g(4<i — x)\ and ip^{x) = a5Cos[k(x — Ad)} + (3^sm[k{x — Ad)}. 
Continuity of ip and its derivative at the five interior boundaries gives ten equations, which 
after much algebra can be reduced to an equation for the energy eigenvalue: 

u D j(£; 4 - k 2 q 2 + q A ) tan 2 (fed) tanh 2 (gd) + fc 2 g 2 (3 + tanh 2 (gd) 

-t&n 2 (kd)) + Akq(q 2 - k 2 ) t&n(kd) tanh(grf) j = 0. (8) 

we again factor out the one period solution giving, 

(A; 4 - k 2 q 2 + g 4 ) tan 2 {kd) tanh 2 \qd) + fc 2 g 2 (3 + tanh 2 (grf) - tan 2 (A;d)) 

+ Akq(q 2 - k 2 ) tan(kd) tanh(qd) = 0. (9) 

As expected from the two period case, Eq. ^ has slightly lower energy eigenvalues than the 
two period case, Eq. (|6]), and there are valid solutions for lower values of V than the two 
period case, as seen in Fig. [2] for the points labeled as three period solution. These solutions 
are lower in energy than the one period Dirichlet solutions, which can be understood by 
looking at the wavefunctions in Fig. |6j The solution for the one period eigenvalue shows 
how the wavefunction is replicated and reduced as x/d increases. The solution for the three 
period equation, Eq. ([9]), puts the majority of the probability in the second and third wells, 
in contrast to the one period solution which has most of the probability in the first well. 
Thus, as in the two period case, the lowest energy eigenvalue is lower, but there remains a 
state whose energy eigenvalue remains unchanged as the number of periods increases. 

For the Neumann boundary condition the wavefunctions for the first and last layer are 
changed to ipi(x) = cticos(kx) and ipQ = a& cosh[g(6<i — x)], with all the others from the 
Dirichlet condition unchanged. The continuity of the wavefunction and its derivative can 
be used to produce ten equations which can be reduced to the following equation for the 
energy eigenvalue, 

u N \{k A ~ k 2 q 2 + q 4 ) tan 2 (kd) tanh 2 (qd) + &V(3 + tanh 2 (grf) - tan 2 (A;<i)) 

+ Akq(q 2 - k 2 ) taa(JW) tanh(gd) } = 0. (10) 

The solution consists of the product of the one period Neumann condition, along with the 
same equation as in Eq. (|9j). As for the two period Neumann solution, the energy states 
are the one period solution (which is the lowest) followed by the same energies as the three 
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period Dirichlet condition. Comparison of the wavefunctions shown in Fig. [7] shows the one 
period solution, with the amplitude decreasing as the power of the number of periods; the 
three period solution has a node and negative values over much of the range of x/d, which 
would be unrealistic for the energy gap function A. 

From the structure of these solutions it is apparent that for a N period solution, the 
Dirichlet and Neumann boundary conditions have an energy eigenvalue that is independent 
of N, which is the lowest energy eigenstate for the Neumann case, although a proof of this 
property has not been shown here. In addition, these results imply there is a common energy 
eigenstate between the two conditions, which is the lowest energy eigenstate for the Dirichlet 
case. 

VI. APPLICATION TO PROXIMITY EFFECT SYSTEMS 

In this study, the idea that the lowest energy eigenvalue should decrease as the system size 
increases is seen to be strongly dependent on the boundary conditions used. For Dirichlet 
conditions the decrease is observed, while for Neumann conditions, the lowest energy eigen- 
value is unchanged as the system size increases. As stated in Ref. [3] the exact proximity 
equations are similar to a one-dimensional Schrodinger's equation with simple potentials and 
Neumann conditions. The one period system faithfully reproduces the behavior of a simple 
bilayer. The results from this study imply that the transition temperature for a N period 
bilayer system (which correlates to the lowest energy eigenvalue) would be independent of 
(In Ref. m one-period, two-period and three-period structures were measured.) However 
the shape of the resulting wavefunction found here is unexpected. If the minimum value of 
A(x) is the energy gap in the system,^ then as the number of periods increases, either the 
gap would become smaller, or the maximum of the self-energy function would increase for 
a fixed gap, neither of which is reasonable, given that the transition temperature remains 
the same. The actual behavior of such a system is still an open question. The simplicity 
of this approach would allow other choices of potential energy well structures to be studied 
in order for students to model various layered films, such as trilayers, or fractal like layer 
arrangements. 
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FIG. 1: Potential wells for bilayers, shown for (a) one period and (b) two periods. The widths of 
the individual layers is the same. 
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FIG. 2: Energy eigenvalues versus potential step height for the one period Dirichlet and Neumann 
boundary condition solutions, Eqs. ^ and Q, respectively; the two period solution with the 
one period solution removed, Eq. ([6]); and the three period solution with the one period solution 
removed, Eq. Q. Notice the Neumann eigenvalues for one period are always below those for the 
Dirichlet solution for one period or the common solutions for multiple periods. 
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FIG. 3: Wavefunctions versus x/d for the one period system for Dirichlet and Neumann boundary 
conditions. Here V = 5. The energy eigenvalues are E = 4.38 for Dirichlet case and 1.12 for the 
Neumann case. 
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FIG. 4: Comparison of wavefunctions for the two period case (Dirichlet condition) with V = 5 for 
E = 4.38 (one period solution) and E = 2.48 given by Eq. [6]). The one period solution puts the 
bulk of the probability in the first layer, while the two period solution puts the bulk in the third 
layer, allowing for a lower energy eigenvalue. 
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FIG. 5: Comparison of wavefunctions for the two period case (Neumann condition) with V = 2 for 
E = 0.70 (one period solution) and E = 1.51 given by Eq. (|6]). Here the two period solution forces 
the wavefunction to be negative over most of x/d. 
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FIG. 6: Comparison of wavefunctions for the three period case (Dirichlet condition) with V = 5 
for E = 4.38 (one period solution, Eq. Q) and E = 2.22 (three period solution, Eq. ([9])). Similar 
to the two period case, the three period solution has the bulk of the probability in the third and 
fifth layers. 
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FIG. 7: Comparison of wavefunctions for the three period case (Neumann condition) with V = 2 
for E = 0.70 (one period solution, Eq. Q) and E = 1.18 (three period solution, Eq. (|9|). As in 
the two period case, the three period solution is negative over most of x/d again. 
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